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Abstract 



> 

^^ , We report a study of T and CP violation in B^B^ mixing using an inclusive dilepton sample 

r — ' collected by the BABAR experiment at the PEP-II B Factory. The asymmetry between £+€+ 

(^ , and i~(.~ allows us to compare the probabilities for B^ -^ B^ and B^ -^ B^ oscillations and 

thus probe T and CP invariance. 

A sample of 20,381 same-sign dilepton events is selected in the 1999-2000 data sample, 
^ ' corresponding to an integrated luminosity of 20.7 fb~ on the T{AS) resonance. We measure a 

D , same-sign dilepton asymmetry of At — (0.5 ± 1.2 ± 1.4) %; for the parameter eb representing 

O ; T and CP violations in mixing, we obtain a preliminary result of 

(D : 

'^ ■ ^<£b) = (1.2 ± 2.9 ± 3.6) X 10-3. 
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1 Introduction 

Since the first discovery of CP violation in 1964 |j^], the kaon system has provided many other 
results probing the CPT and T discrete symmetries Q- The BABAR experiment is not limited to 
an investigation of CP violation through the measurement of sin(2/?) or the angle a; in a similar 
way as for kaon system studies, it is also possible to investigate CP violation purely in mixing and 
disentangle whether this CP violation is due to T or CPT violation. 

In this document, we have adopted a formalism very similar to that used for the kaon system Q. 
In the absence of a CP phase, we can define the CP eigenstates \Bi) and I-B2) ^^ 

\Bl) = -^(|S°) + |5°)), 

\Bl) = ^^{\B')-\W)). 

The physical states (solutions of the effective Hamiltonian) can be written as 

1^0) = , ,^ ^ [\Bl) + {eB + 6B)\Bl)l 
^J\ + \eB + SbV 

\B%) = ^ j B^^) + ^eB-5B)\B',)]. 

In the case of CPT and CP invariance, 6b is equal to 0. Similarly, T and CP invariance gives 
eb = 0. This means that CP violation in mixing requires either eb "/" or 6b j^ 0. Inclusive 
dilepton events in BABAR provide a very large sample with which to study CP violation in mixing 
and test T and CPT conservation. The semileptonic (muon or electron) branching fraction of B 
mesons is about 20%. Therefore, dilepton events represent 4% of all T(45) -^ BB decays. The 
flavor of the B is tagged by the sign of the lepton. Assuming CP invariance in direct semileptonic 
decays, the asymmetry between same-sign dilepton pairs, i~^i~^ and i~i~ , allows a comparison of 
the two oscillation probabilities P{B^ — > B^) and P[B^ — > B^) and therefore probes T and CP 
invariance: 

^ (^^. ^ P{B^ ^ B^,At) - P{B^ ^ B^,At) ^ 4Re(£B) ,^. 

^^ ' P(BO^BO,At) + P(BO^BO,At) ^ l + ksP' ^^ 

This asymmetry^ does not contain the CPT violation parameter 6b to first order. For this asymme- 
try to be different from zero both T and CP violation in mixing are required. In the approximation 
that |ri2| << IM12I, At = Im(ri2/Mi2) where M12 — 2 ^^12 is the off-diagonal element of the 
complex effective Hamiltonian for the B^- B^ system. Standard Model calculations [^ predict the 
size of this asymmetry to be of order (0.5 — 5.0) x 10"'^. Therefore, a large measured value for this 
asymmetry could be an indication of new physics. 

The measurement of Aj- reported here is performed using events collected by the BABAR detector 
at the PEP-II asymmetric B Factory between October 1999 and October 2000. The integrated 
luminosity of this sample is 20.7 fb~ taken on the T{4S) resonance ( "on-resonance" ) and 2.6 fb~ 
taken 40 MeV below the resonance ("off-resonance"). The BABAR detector and its performance 
are described elsewhere Q. 

^In another formalism M], the physical states are related to the S" flavor eigenstates by |-B£) = p\B'^) + q\B ) 
and \B%) = pI-B") — g|S") where p and q are complex mixing parameters with the normalization |pp + |g[^ = f . The 
charge asymmetry in terms of \q/p\ is therefore equal to At = (1 ~ \q/p\'^)/{l + \q/p\^). 



The organization of this paper is as follows. The particle identification and event selection are 
described in Sec. |2|. In particular, this section shows how cascade leptons from charm decays from 
charged B or unmixed neutral B events are rejected with a neural network (NN) approach and a 
Az cut at 200 ^m. The method to correct the charge asymmetry in the detection of the leptons 
is explained in Sec. ^. Finally, Sec. |^ shows the details of the fit on data and gives the evaluation 
of systematic uncertainties. To avoid the possibility of bias, both the time distribution of the 
charge asymmetry and the number of positive and negative same-sign dileptons are blinded during 
the analysis. The unblinding of the At measurement is performed once all studies of systematic 
uncertainties are finished. 

2 Selection of dilepton events 

In this study of T and CP asymmetries, the flavor of the B meson at the time of its decay is 
determined by the sign of direct leptons produced in semileptonic B decays. This section describes 
the selection of leptonic tracks and the rejection of cascade leptons. 

2.1 Lepton identification 

Lepton candidates must have a distance of closest approach to the nominal beam position in the 
transverse plane less than 1 cm, a distance of closest approach along the beam direction less than 
6 cm, at least 12 hits in the Drift Chamber (DCH), at least one z-coordinate hit in the Silicon 
Vertex Tracker (SVT), and a momentum in the T(4S) center-of-mass system (CMS) between 0.7 
and 2.3 GeV/c. 

Electrons are selected by specific requirements on the ratio of the energy deposited in the 
Electromagnetic Calorimeter (EMC) and the momentum measured in the DCH, on the lateral 
shape of the energy deposition in the calorimeter, and on the specific ionization density measured 
in the DCH. Muons are identified through the energy released in the calorimeter, as well as the 
strip multiplicity, track continuity and penetration depth in the Instrumented Flux Return (IFR). 
Lepton candidates are rejected if they are consistent with a kaon or proton hypothesis according 
to the Cherenkov angle measured in the Detector of Internally Refiected Cherenkov Light (DIRC) 
and to the ionization density measured in the DCH. 

The performance and charge asymmetry of the lepton selection are determined with data con- 
trol samples, as a function of the particle momentum as well as the polar and azimuthal angles. 
The electron and muon selection efficiencies are about 92% and 75%, with pion misidentification 
probabilities around 0.2% and 3%, respectively. All corrections of the charge asymmetry in the 
identification of leptons are discussed in Sec. §. 

2.2 Background rejection 

Non-i?i? events are suppressed by requiring the Fox- Wolfram ratio of second to zeroth order mo- 
ments H to be less than 0.4. In addition, the residual contamination from radiative Bhabha and 
two-photon events is reduced by requiring the squared invariant mass of the event to be greater 
than 20 ( GeV/c^)^, the event aplanarity to be greater than 0.01, and the number of charged tracks 
to be greater than four. 

Electrons from photon conversions are identified and rejected with a negligible loss of efficiency 
for signal events. Leptons from J/tp and ipC^S) decays are identified by pairing them with other 
oppositely-charged candidates of the same lepton species, selected with looser criteria. We reject 



the whole event if any combination has an invariant mass within 3.037 < M{i'^£ ) < 3.137 GeV/c^ 
or 3.646 < M{£+£-) < 3.726 GeV/c^. 

2.3 Selection of direct dileptons 

To minimize the dilution due to wrong-sign leptons from cascade charm decays (produced in 6 ^ 
c ^ I transitions), we use a NN algorithm that combines five discriminating variables. These are 
calculated in the CMS and are 

• the momenta of the two leptons with highest momentum, p| and j>2i 

• the total visible energy Etot and the missing momentum Pmiss of the event; and 

• the opening angle between the leptons, 9i2- 

The distributions of these variables are shown in Fig. |l], for data and Monte Carlo simulation. 
The first two variables, pi and P2, are very powerful in discriminating between direct and cascade 
leptons. The last variable, 9i2, efficiently removes direct-cascade lepton pairs coming from the same 
B and further rejects photon conversions. Some additional discriminating power is also provided 
by the other two variables. The NN architecture (5:5:2) consists of three layers, with two outputs 
in the last layer, one for each lepton in the event. The network is trained with 40,000 dileptons 
from generic Monte Carlo B^B^ and B^B^ events. The outputs are chosen to be 1 and for direct 
and cascade leptons, respectively. The same network is used for both electron and muon selection. 
We require both outputs to be greater than 0.8. 

2.4 Background rejection based on Az information 

In the inclusive approach used here, the z coordinate of the B decay point is the z position of 
the point of closest approach between the lepton candidate and an estimate of the T(45) decay 
point in the transverse plane. The T{4S) decay point is obtained by combining the beam spot 
constraint and the relative position of the two lepton tracks. The time difference At between the 
two B mesons is determined from the difference in z between the two B^s by At = Az/{P^)c, where 
ip-f) w 0.56. 

To measure At, we want to select mixed events {B^B^ or B^B^ pairs) whose time depen- 
dence varies as e~ '"^bo [i — cos(AmrfAt)]. The main sources of background (cascade leptons from 
unmixed B^B^ events and B^B~ events) vary respectively as ~ e~ ''^bo [i -\- cos(AmrfAt)] and 
^ g-At/Tg± (see Fig. |^. Therefore, a requirement of Az > 200 ^m allows us to eliminate a large 
fraction of background without dramatically decreasing the signal efficiency. A Az cut is also effec- 
tive at removing backgrounds such as non-BB events or J/tp decays. Finally, in the measurement 
of At, the dilution factor due to remaining background will be corrected as a function of Az. 

2.5 Event yields and sample composition 

Application of the selection criteria described above results in a sample of 20,381 same-sign dilepton 
events, consisting of 5,252 electrons pairs, 5,152 muons pairs and 9,977 electron-muon pairs. 

The fraction of non-BB events, measured with the off-resonance data, is 4.3% with a charge 
asymmetry of (—5 it 10)%. 
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Figure 1: Distributions of the discriminating variables (a) pi, (b) P2, (c) Etot, (d) Pmiss and (e) 6*12, 
for data (points) and Monte Carlo events (histograms). The contributions from direct-direct pairs, 
direct-cascade or cascade-cascade pairs, and pairs with one or more fake leptons are shown for the 
Monte Carlo samples. 



11 



0.006 



0.005 



0.004 



0.003 



0.002 - 



0.001 




200 400 600 800 1000 1200 1400 1600 1800 2000 

Az (^m) 



Figure 2: Probability density functions determined from a fit to the data with the value of Arrid fixed 
to the world average value [^ , for the signal [B^B^ or B^B^ pairs) (solid line) and the background 
(cascade from unmixed B^B^ and B^B" events, and non-BB events) (dotted line) as a function of 

Az. 

3 Detector-induced asymmetries 

Since the asymmetry At is expected to be small, we must control possible charge asymmetries 
induced by the detection and reconstruction of electrons and muons. The strategy in this analysis 
is first to correct on an event-by-event basis charge asymmetries determined with independent 
control samples, and then to validate this approach with a sample of single direct leptons from B 
decays. 

3.1 Charge asymmetry correction 

The three sources of charge asymmetry in the selection of lepton candidates are 

1. the difference in tracking efficiency for positive and negative particles, e^^^^, 7^ ^track'^ 

2. the difference in particle identification efficiency for positive and negative leptons, ej"^^ / e~^^\ 
and 

3. the difference in misidentification probability for positive and negative particles, •q'^^^ 7^ 'Hnid- 

These efficiencies and probabilities are estimated using independent samples (see next sections) on 
an event-by-event basis and as a function of several kinematic variables Xi of the charged track: the 
total momentum, transverse momentum, polar angle and the azimuthal angle of the charged track. 
The numbers of "detected" positive and negative leptons ^detected^^) ^'^^ related to the numbers of 
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true leptons N^^^^{i) by the equation: 

+r{K,p*) ■ r]^.^{K,Xi) + r{p,p*) ■ r]^.^{p,Xi)], 

where r{TT,p*), r{K,p*) and r(p,p*) are the relative abundances of hadrons (vr, K, and p) with 
respect to the lepton abundance for a given p* (the momentum of the track in the CMS). These 
quantities are obtained from generic BB Monte Carlo events, after applying the event selection 
criteria. To correct for the charge asymmetries, we apply a weight proportional to the ratio 
^truei^) / ^detected(^) ' ^^^ ^^^^ lepton in the sample. 

3.2 Correction with control samples 

3.2.1 Charge asymmetry in tracking 

The event selection requires at least 12 DCH hits per track, which can introduce a charge asymmetry 
in the tracking efficiency. In this analysis, we are selecting tracks with large momentum in the CMS, 
which implies that the transverse momentum is also large. Therefore, the dilepton sample should 
be only slightly affected by any charge asymmetry in tracking. To remove this possible bias, we 
determine separate efficiencies for positive and negative particles. The tracking efficiency, which 
is dominated by the DCH, is defined as the ratio of the number of SVT tracks with 12 DCH hits 
divided by the initial number of SVT tracks. These track efficiency correction tables are computed 
as a function of transverse momentum, polar angle and azimuthal angle in the laboratory frame. 
The charge asymmetry correction is less than 0.1% on average in the relevant ranges for the lepton 
tracks. 

3.2.2 Charge asymmetry in lepton identification 

Since the particle identification efficiencies and the misidentification probabilities are determined 
as the ratio of the number of events with an identified track over the number of initial events, the 
efficiency and misidentification probability tables are independent of the tracking efficiency. These 
tables are computed as a function of total momentum, polar angle and azimuthal angle in the 
laboratory frame for different control samples: 

• The identification efficiencies for the electrons are measured with the combination of two 
control samples: 77 -^ eeee and radiative Bhabha events; 

• The identification efficiencies for the muons are measured with a control sample consisting of 
77 -^ ee/U/i events; 

• The misidentification probabilities are determined using control samples of kaons produced 
in D*^ -^ 7r~^D^ -^ tt'^K^tt'^ decays (and charge conjugate), pions produced in Ks -^ vr+vr" 
decays, one-prong and three-prong r decays, and protons produced in A decays. 

For the electrons, the charge asymmetry in the efficiency reaches (0.5-1.0)% in some regions of the 
lepton phase space. The impact of the charge asymmetry in misidentification is negligible because 
the absolute misidentification probability for pions is extremely small (~ 0.2%). However, the A 
control sample indicates a very large misidentification probability for antiprotons with momentum 
~ 1 GeV/c. Such an effect is due to the annihilation of antiprotons with protons in the calorimeter, 
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which produces a signature similar to that of an electron. The impact of this effect is balanced by 
the low relative abundance of antiprotons in generic B events. Overall, antiprotons induce a charge 
asymmetry of order 0.1% and a correction is applied for this effect. 

For the muons, the ee///i control sample shows that the charge asymmetry in the efficiency 
reaches 0.5%. The fraction of fake pions (~ 3%) is much larger than in the case of electrons but 
there is no indication of any charge asymmetry. On the other hand, the kaon misidentification 
distribution shows a charge asymmetry at the level of (10-20)% due to the difference between 
cross sections for K^ and K~ mesons interactions with matter in the range of momentum around 
1 GeV/c. 

3.3 Validation with single direct leptons 

A cross check of the correction for charge asymmetries in the lepton selection is performed with 
an independent sample that has the same topology and kinematics as the leptons from dilepton 
events: namely direct leptons from semileptonic B decays. The selection of these events is quite 
similar to the dilepton selection described in Sec. |2[ 

The single-lepton charge asymmetry is sensitive to the charge asymmetry due to detection bias 
but it may also be affected by the real physical charge asymmetry At in the dilepton events. In 
this case, the physical charge asymmetry A^i^aie^ '"^'^ ^^ written as 

where D is a dilution factor due to background (cascade leptons, etc.), Xd is the B^ mixing param- 
eter, and R is the fraction of charged B mesons in the sample.^ The possible bias introduced by 
At is suppressed by more than one order of magnitude and is therefore neglected. 

With the 1999-2000 data set, we select roughly 1.5 million electrons and 1.5 million muons. 
After subtraction of scaled off-resonance data and after applying a correction weight as defined in 
Eq. H, we measure the remaining asymmetry as a function of the total momentum (see Fig. |3|). For 
both muons and electrons, the distributions are consistent with being flat. The single-muon charge 
asymmetry (— 0.35ib0.17)% and the single electron charge asymmetry (— 0.30ib0.14)% are consistent 
with zero. These measured charge asymmetries for single leptons show that the systematic errors 
related to the muon and electron detection are at the level of ±0.35% and ±0.30%, respectively. 

4 Measurement of Re(£B) 
4.1 Measurement of At 

Equation |I| is applicable for pure signal [B^B^ and B^B^ pairs). However, the dilepton events 
are contaminated by cascade leptons from B^B~ and unmixed B^B^ events (see Fig. |2|), non-BB 
events, and J/ij: decays. Assuming no charge asymmetry in the background and assuming CP 
invariance holds in the direct leptonic B decays,^ we can write the measured asymmetry A^^""^ 

^ R is the ratio (6+/h )/(6o/oo) where 6+ and bo are respectively the semileptonic branching fractions of charged 

and neutral B mesons, and /+_//oo is the production ratio of charged and neutral B mesons. 

^In the literature, the equality of the decay probabilities PIB*^ — > £'^) and P{B^' — > £^) is usually obtained 
from CPT invariance in the decay for inclusive production of direct leptons. However, in this analysis, the signal is 
selected by imposing cuts on the lepton momentum, which is equivalent to considering a partial decay channel. As 
a consequence, the equality of the partially integrated spectrum for leptons and antileptons requires CP symmetry, 
which is a stronger assumption. 
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Figure 3: Single direct lepton control sample: charge asymmetries for electrons (left plot) and 
muons (right plot) as a function of the total momentum. 



(see Fig. ^ in terms of the number of events A^ as 



XAt) 
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A' 



T ■ 



S{At) 



S{At) + B{At)' 



(3) 



where S{At) is the number of signal events and B(At) the total number of background events. 
The assumption of no charge asymmetry in the background is confirmed by the off-resonance data 
where the charge asymmetry (—5 it 10)% is consistent with zero. In addition, the charge asymmetry 
of the events with Az < 100 /xm, which contain 85% background (cascade leptons from B^ and 
unmixed B^), is (1.2 it 1.4)% and so is also consistent with zero. Finally, a possible dilution of Aj^ 
due to double mistag is neglected because the probability of double mistag is at the level of only 
1%. All these assumptions are taken into account in the determination of the systematic error. 

Therefore, extraction of a value for At requires a determination of the dilution factor 
S{At)/[S{At) + B{At)]. The dilution factor can be measured directly from the data with the 
probability density functions (Fig. g) obtained from the full dilepton sample with the value of 
Anid fixed to the world average value Q. With this method, the fraction of non-BB events is 
determined from off-resonance data, and the fraction of cascade leptons and the resolution function 
corrections are measured directly from the dilepton data. In order to benefit in the determination 
of At from the full information provided by Az, a correction factor for the background dilution 
1+ B (Azi) / S (Azi) is applied to each Azi bin considered in the fit. 

After applying the above correction to the dilepton sample, we measure At = (0.5 it 1.2) % 
from a fit to the distribution of the charge asymmetry as a function of At (see Fig. Q) . 
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Figure 4: Corrected dilepton charge asymmetry as a function of At. The hne shows the result of 
the fit for the same-sign dileptons with Az > 200 //m. 

4.2 Systematic uncertainties 

The detection charge asymmetry is partially corrected by applying an event-by-event weight. We 



assign the residual asymmetries measured with the single lepton samples (see Sec. 3.3), ±0.30% 
for the electrons and ±0.35% for the muons, as the systematic errors due to charge asymmetry 
in detection efficiencies. After taking into account the dilution factor, the total systematic error 
related to the charge asymmetry in the detection efficiencies are ±0.5% and ±0.6% for electrons 
and muons, respectively. 

The charge asymmetry of (—5 ± 10)% measured with the off-resonance data leads to a ±0.7% 
uncertainty on the At measurement, determined from the statistical error ±10%. In a similar way, 
the charge asymmetry of (1.2 ± 1.4)%, obtained with the events satisfying Az < 100 ;um, results 
in a ±0.9% uncertainty on At- If we assumed CP invariance in the decays producing the cascade 
leptons, this uncertainty would vanish. 

The other systematic uncertainties are due to the background dilution correction. This correc- 
tion is measured with the data from the full dilepton sample with the value of Am^ fixed to the 
world average value |p. The uncertainty on the ratio B/S leads to a ±3% multiplicative error on 
At, which is negligible. 

All sources of systematic uncertainty are listed in Table ||. The total systematic uncertainty is 
±1.4%. 



5 Conclusions 

With the 1999-2000 data consisting of 20.7 fb^ on-resonance and 2.6 fb^ off-resonance, we have 
selected 20,381 same-sign dilepton events with a NN approach and a rejection of the cascade leptons 
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Type of systematic error 


ct{At){%) 


Electron charge asymmetry in the detection 
Muon charge asymmetry in the detection 
Non-BB background charge asymmetry 
BB background charge asymmetry 


0.5 
0.6 
0.7 
0.9 


Total 


1.4 



Table 1: Summary of systematic uncertainties on At- 



from charged B or unmixed neutral B events based on the Az information. Charge asymmetries in 
the lepton selection are corrected event-by-event and the correction method is confirmed with the 
single lepton sample. We have measured a same-sign dilepton asymmetry of A^ = (0.5ibl.2ibl.4)%. 
From Eq. ||, the At asymmetry gives a preliminary result for the T and CP violation parameter 

— ^, = (1.2 ± 2.9 ± 3.6) X 10"^. 

This preliminary measurement is the most stringent test of T and CP violation in B^B^ mixing to 
date and is consistent with previous measurements [^ of Re(eB)/(l + |£Bp)- With the formalism Q 
using the complex parameters p and q, this measurement of At gives \q/p\ = 0.998 ± 0.006 ± 0.007. 
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